arXiv: 1506.0633Ivl [cond-mat.other] 21Jun2015 


Motional Averaging of Nuclear Resonance in a Field Gradient 

Nanette N. Jarenwattananon and Louis-S. Bouchard* 

Department of Chemistry and Biochemistry, University of California Los Angeles, 

607 Charles E Young Drive East, Los Angeles, CA 90095-1059 
(Dated: June 23, 2015) 

The traditional view of nuclear-spin decoherence in a field gradient due to molecular self-diffusion 
is challenged on the basis of temperature dependence of the linewidth, which demonstrates different 
behaviors between liquids and gases. The conventional theory predicts that in a fluid, linewidth 
should increase with temperature; however, in gases we observed the opposite behavior. This 
surprising behavior can be explained using a more detailed theoretical description of the dephasing 
function that accounts for position autocorrelation effects. 

PACS numbers: 05.20.Jj, 51.10.-|-y, 51.20.-|-d, 66.10.cg, 82.56.Lz 


For over six decades, diffusion-weighted nuclear mag¬ 
netic resonance (NMR) has been the flagship experi¬ 
ment for measurements of molecular self-diffusion in free 
or confined geometries. Diffusion-based NMR experi¬ 
ments have a wide range of applications from porous me¬ 
dia [1], catalysis [2], materials science and chemistry [3] to 
biomedicine [4]. Consider the simple experiment shown 
in Fig. 1, where the nuclear induction signal is read out 
in the presence of a magnetic-held gradient. The gra¬ 
dient modulates the magnetization spatially along the 
gradient’s direction (assuming a sufficiently strong exter¬ 
nal held so that the gradient is unidirectional). Time 
evolution of this magnetization in the presence of dif¬ 
fusion effects provides a direct and unambiguous mea¬ 
surement of the self-diffusion process. In the traditional 
description of molecular self-diffusion [5-9], a molecule 
undergoes a random walk whereby at each time step, 
the nuclear spins accumulate phase increments that are 
randomly drawn from a normal distribution. In the pres¬ 
ence of a magnetic-held gradient, the decoherence of the 
nuclear induction signal S{t) follows the well-known text¬ 
book expression [5-10]: 

S{t) = exp (1) 

where 7 „ is the nuclear gyromagnetic ratio, D is the self¬ 
diffusion coefficient, g is the applied gradient strength, 
and t is time. The dependence has been extensively 
validated and has been utilized to measure molecular self¬ 
diffusion coefficients in a wide variety of liquids [10-12]. 
In a gas, however, the situation is more complicated. The 
assumption of a normally-distributed phase accumula¬ 
tion at every time step is difficult to justify in light of the 
fact that gas molecules undergo much more rapid motion 
than in liquids, due to much longer free displacements be¬ 
tween collisions. The farther the molecular displacements 
along the direction of the gradient, the faster the nuclear 
spins will lose memory of their immediate environment. 
This memory loss should be expected to enter the de¬ 
scription of the decoherence process. Surprisingly, this 
simple aspect of free diffusion is still lacking a thorough 
experimental verification. 
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FIG. 1. Measurement of molecular self-diffusion in a constant 
field gradient. A 90° radiofrequency (rf) pulse tips the magne¬ 
tization. This resulting nuclear induction signal is measured 
in a constant gradient of amplitude g. The time evolution of 
magnetization, as described by the textbook [10-12] expres¬ 
sion (1), interrogates the molecular self-diffusion process. 


The NMR signal, S{t), from an ensemble of spins ini¬ 
tially located at x(0) is given by the expectation value of 
the phase factor: 


S{t) = ^exp J , (2) 


where uj(t) describes the time-dependent resonance fre¬ 
quency offset (in the rotating frame). In the presence of 
a magnetic field gradient g, the resonance frequency is 
uj{t) = ^ngxit), where x{t), the position of the nuclear 
spin after time t, is a random process. If we assume a 
Gaussian random process that is stationary in the wide 
sense, this expectation value takes the form [9]: 

exp J {x{t'))dt' - j^g"^ J {x{t')x(0)){t - t')dt' 

( 3 ) 

The first term, exp(* 7„5 Ji^{x{t'))dt'), encodes the posi¬ 
tion in the phase of the spins [13] after undergoing dis¬ 
placement but does not lead to decoherence unless the 
diffusion is not free. The second term describes signal 
decay, and consequently, affects linewidth. More compli¬ 
cated random processes may lead to higher-order terms, 
but we shall limit our discussion to the case where x{t) 
is a wide-sense stationary Gaussian random process. 
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Over the timescale t of the NMR measurement, liq¬ 
uid phase molecules experience displacements that are 
much smaller than those in the gas phase. We may write 
this as x{t) « a;(0), for a liquid. This approximation, 
known as the Einstein-Fick limit, indicates that the po¬ 
sition autocorrelation function can be approximated by 
the mean-square displacement: 

{x{t)x{0)) « {x{t)x{t)) = 2Dt, (4) 

where D is the self-diffusion coefficient. It is known 
from experiments that in liquids, this limit holds [10- 
12]. However, Duquid and Dgas differ by three orders 
of magnitude, and it is unclear if this approximation 
also holds for gases. Away from the Einstein-Fick limit, 
contributions from ballistic transport become significant. 
We note that substitution of the limit (4) into the second 
term of (3) recovers the result (1) as a special case. 

In this Letter, we show that (1) does not hold for gases 
based on the analysis of linewidth as function of temper¬ 
ature (T). Because D{T) increases with T, (1) predicts 
that linewidth should increase with T. We find that gases 
instead undergo line narrowing with temperature. It is 
unclear by cursory inspection of (1) how gases may differ 
from liquids, especially in light of the fact that the de¬ 
pendence of line broadening on diffusion coefficient has 
been verified experimentally in several studies (see, for 
example, [11, 14, 15]). The key to establishing this dis¬ 
tinction is a closer look at the temperature dependence 
of the line broadening mechanism, as explained below. 

In the NMR experiment, nuclear spins are well isolated 
from the lattice and do not depolarize or randomize their 
phases when undergoing molecular collisions, in contrast 
to collisional broadening mechanisms in optics. Thus, the 
description of line broadening in such a “weak collision” 
regime reflects the histories of molecular displacements. 
The simplest way to account for this is through a posi¬ 
tion autocorrelation function. Suppose that the particle 
displacements are modeled using a generalized Langevin 
equation (GLE) with memory kernel: 

Mv + f T{t — t')v{t')dt'= r]f{t), (5) 

Jo 

where M is the mass of the diffusing particle, r(t) is a 
memory kernel, v{t) = x{t) is the particle velocity, v is 
its acceleration, and ? 7 /(t) is a stochastic force. The GLE 


has been validated experimentally for Brownian particles 
(M 3> TO, where to is the mass of fluid particles); for ex¬ 
ample, in the studies [16-19] M was 10^° times larger 
than TO. So while the GLE was not designed to model 
self-diffusion processes, it can be invoked to model vis¬ 
cous drag effects via the memory kernel. In what follows, 
we shall set M = 10^*^ - to, which is the only regime we are 
aware of, where the GLE has been validated experimen¬ 
tally based on direct measurements of individual histories 
(namely, in Refs. [16-19]). 

By the fluctuation-dissipation theorem, the stochas¬ 
tic force ??/(t) describes colored noise, (? 7 /( 0 )? 7 /(t)) = 

kTV{t), where k is Boltzmann’s constant. The 

time-correlation function {x{t)x{0)) is obtained from 
{v{t)v{0)) by integrating twice the velocity autocorrela¬ 
tion function 

{v{t)v{0)) = -^{x{t)x{0)). (6) 

Projecting equation (5) with the operator (v(0), •) yields 
the deterministic equation: 

M{v{0)v{t)) + f T{t — t'){v{0)v{t'))dt' = 0. ( 7 ) 

Jo 

Integrating this velocity autocorrelation function once 

Ht) = f\v{0)v{t'))dt', ( 8 ) 

Jo 

and using the equipartition theorem, (v(O)u(O)) = 

kT/M^ as the initial condition, we get: 

Mv{t) + f T{t — t')v{t')dt' = kT. (9) 

Jo 

For the memory kernel to describe the delayed response of 
the surrounding fluid, we choose the Ornstein-Uhlenbeck 
process: 

r(t) = ( 7 ^/to) exp(—yt/m), 

where 7 is the friction coefficient which is proportional 
to the viscosity of the medium and to is the mass of 
molecules in the surrounding medium causing friction. 
Denoting = ( 7 / 2 to )(1 =F \/l — Am/M), we obtain 
the solution to equation (9): 
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( 10 ) 


From this we get the position autocorrelation function. 
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(a;(t)a;(0)) 


kT 

M(C+-C-) 



This result was also derived by N 0 rrelykke [20] using a 
different method. In the Einstein-Fick approximation, 
a;( 0 ) « x(t), and this position autocorrelation function 
reduces to 2Dt. Equation (11) generalizes S{t) outside 
the Einstein-Fick limit. There are three distinct regimes: 
overdamped (M > 4m), critically damped (4m = M), 
and underdamped (M < 4m). Standard Brownian mo¬ 
tion of large particles is strongly overdamped (M » m). 
The Einstein-Fick limit occurs when the ratio yt/m 
is sufficiently large to cause appreciable decay of the 
Ornstein-Uhlenbeck kernel. Using 7 based on the Stokes’ 
law and t = 40 /iS as the sampling time of the nuclear 
induction signal, we find that typical values of this ra¬ 
tio for liquids are "ft/m ~ 1 whereas for gases we have 
"ft/m < 1 . 

In the case of a gas, we may obtain the overall tem¬ 
perature by modeling 7 using the Stokes’ law, 7 = STrry^d 
{rfy, shear viscosity of the medium; d, sphere diameter), 
which holds in the limit of low Reynolds numbers and 
invoking the Sutherland’s formula [21] for rjy: 


Mo(ro + c)(r/ro)3/2 ^ 

T + C ^T + C" 


( 12 ) 


where C is Sutherland’s constant for the gas, fio is 
the viscosity at temperature Tq. At low temperatures, 
rjy ~ whereas at high temperatures, rjy ~ y/T. By 
substituting this into ( 12 ), we obtain at an expression for 
the envelope function of the signal decay for a gas which 
does not rely on the Einstein-Fick limit: 


S{t) = exp{-"flg^Kt) (13) 


with 
niT) = 


kT{-mCC+ - mC-C+ + C-7 + C-C-h7 + C+7) 


The linewidth A/ follows the power law: 


A/ 


t<c 

T-1/2, T>C 


(14) 

(15) 


which predicts a temperature dependence that is opposite 
(i.e., line narrowing with increasing temperature) to that 
based on self-diffusion in the Einstein-Fick limit (1). An 
analogous expression in the case of liquids can be derived 
using a suitable model for the temperature dependence 
of the viscosity in a liquid. However, this will not be 
needed here, because we shall see that the linewidth is 
essentially independent of temperature. 


mC+ / 
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-C+i 




( 11 ) 


Measurements of A/ were carried out as a function 
of T for three gases in the high temperature regime 
(T > C, as determined by the Sutherland’s constant for 
each gas [22-25]). The results are shown in Figure 2. 
The average exponent was found to be —0.47 ± 0.04, in 
agreement with the theoretically predicted value of — 1/2 
in Eq. (15). The low temperature regime {T < C) could 
not be investigated due to experimental limitations of 
our instrument. A temperature dependence of linewidth 
could not be detected within experimental error for liq¬ 
uids, as shown in Figure 3, where we investigated nine 
different liquids over the range 180-450 K. We note that 
the Sutherland’s formula (and therefore. Equation 13) is 
applicable to gases only, so a lower limit on temperature 
for the liquids is imposed by the freezing points. 

18 T 
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FIG. 2. Validation of the high-temperature 

regime {T > C) for gases. The values of luT shown corre¬ 
spond to the temperature range T = 180 — 490 K. Three 
different gases were investigated: methane, acetylene, and 
propylene. The temperature dependence on linewidth (scaled 
to gradient strength) was found to be A/ oc j’-0'4':’±o.04 
eraged over the three gases). Although different values of 
the applied gradient g are shown here to avoid overlapping of 
the curves (methane, g=Q.l^ G/cm; acetylene, <7=0.07 G/cm; 
propylene, <7=0.1 G/cm), the scaled linewidth is independent 
of g. 


We now turn out attention to the gradient dependence 
of the line broadening, which, according to (13), should 
be proportional to . In experiments, however, we found 
two regimes: in the limit of weak gradients. A/ oc 
whereas for strong gradients. A/ oc g^ (see Figs. 4 and 
5). The < 7 ^ gradient dependence is predicted according to 
(13); however, the A/ ex < 7 ^ is not. Moreover, these two 
regimes apply to both liquids and gases, according to the 
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FIG. 3. Temperature dependence of linewidth (A/) in liq¬ 
uids. Nine different liquids were investigated, as shown by 
the different symbols. The valnes of ln(T) shown span the 
temperatnre range T = 180 — 450 K. A) For a hxed gra¬ 
dient strength of 3=0.05 G/cm, all linewidths were broad¬ 
ened by a similar amonnt, hence the overlap in the data. 
At fixed g, the linewidth did not exhibit any detectable de¬ 
pendence on temperatnre. B) Increasing applied gradient 
strength did not alter the independence of linewidth on tem¬ 
peratnre. Applied gradient strengths were: nitromethane 
(1 G/cm), dichloromethane (dcm, 0.5 G/cm), acetonitrile 
(0.4 G/cm), chloroform (0.3 G/cm), benzene (0.3 G/cm), wa¬ 
ter (0.2 G/cm), trifluoroacetic acid (tfa, 0.1 G/cm), dimethyl 
sulfoxide (dmso, 0.1 G/cm), acetone (0.05 G/cm). Symbol in 
A refers to the same liquid as in B. 


results of Figs. 4 and 5. The emergence of the regime 
is likely due to the convolution of the line shape with 
the sample shape that arises under an applied gradient 
and forms the basis for frequency encoding in magnetic 
resonance imaging. The full signal equation integrated 
over the sample shape is: 

S{t) = j exp{ijngxt) exp{—ag^t)p{x)dx, (16) 

where p{x) describes the spin density profile along the 
X direction (sample shape), ex.p(i^ngxt) is the frequency 
encoding and exp(—a^^t) is the line broadening accord¬ 
ing to (13). For the particular case where the sample 
shape p{x) is Gaussian (a reasonable approximation for 
the sensitivity profile of a saddle coil, such as the one 
used in these experiments), the Fourier transform of (16) 


FOV (cm) 


1057 236 87 32 12 3 0.6 
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FIG. 4. Dependence of linewidth (A/) on applied gradient 
strength (g) for gases. Three different gases were investi¬ 
gated. Two different regimes are found: in the limit of strong 
applied gradients. A/ scales as whereas for weak gra¬ 
dients A/ scales as All data was acquired at ambient 

temperature. 


with respect to time is a convolution: 

S{u}) = Gauss G Lorentz = Voigt, (17) 

in which the Gaussian is ~ exp(—x^/ 2 ( 7 ^), and the 
Lorentzian is ~ r^/(a:^ -h F^). In terms of the fnll width 
at half maximum of the Gaussian {fc = 2ij'fgy/2 In 2) 
and Lorentzian {Jl = 2r = 2ag^) profiles, the width 
of the Voigt profile can be expressed as /„ « 0.5 /l + 
^0.2/1 + /^. Thus, in the frequency encoding regime, 
where sample shape effects dominate, the line broadening 
behaves as g^ regardless of whether (1) or (13) are used to 
describe diffusion effects. The frequency encoding regime 
is reached when the field of view FOV=/s/ 7„3 ifs, sam¬ 
pling rate; g, gradient amplitude), becomes comparable 
to the size of the rf-sensitive region (~ 1 cm in our ex¬ 
periments). Depending on the applied gradient strength, 
the experimental FOV ranges from 0.6 to 1460 cm. The 
FOV values corresponding to applied gradient ( 3 ) are in¬ 
dicated in the upper horizontal axes of Figs. 4 and 5, 
where the two regimes, g^ and 3 ^, are indicated. 

We have presented a revised expression for line broad¬ 
ening (13) that not only takes into account the autocor¬ 
relation effects in the diffusion process, but also suggests 
that self-diffusion processes in the NMR experiment may 
be described using a stochastic GLE; at least, as far as its 
temperature dependence is concerned. The GLE (13) en¬ 
ables a convenient description of the memory effects aris¬ 
ing from the viscous drag effects, which were essentially 
missing from the traditional description (1). Such drag 
effects yielded the correct temperature dependence for 
gases and have been used in a recent publication to non- 
invasively map temperatures of gases during catalytic re¬ 
actions [26]. The method could also be useful in the 
validation of heat-transfer models for gas-phase thermal 



5 


FOV (cm) 

1875 690 254 93 34 12 2 



-7 -6 -5 -4 -3 -2-10 1 

in(g) 

FIG. 5. Dependence of linewidth (A/) on gradient strength 
{g) for liquids. Nine different liquids were investigated. Two 
different regimes are found: in the limit of strong applied gra¬ 
dients, A/ scales as i whereas for weak gradients A/ 

scales as All data was acquired at ambient temper¬ 

ature. 

exchange systems, which currently rely on numerical re¬ 
sults from computational fluid dynamics models. Finally, 
we note that since the decay function (13) involves the 
first power of time instead of its third power (1), (13) 
could have implications for the design of dynamic decou¬ 
pling schemes for coherent quantum control [27]. 

See Supplemental Material [22] at [URL to be inserted 
by publisher] for experimental methods, data analysis de¬ 
tails, and sample data. 
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SUPPLEMENTARY MATERIAL (SM) 
SM: EXPERIMENTAL SET-UP 


Sample Preparation 

Methane gas (> 99% purity) was purchased from Air- 
gas, Inc. and used as provided. Acetylene(> 98% purity) 
and propylene gases (> 95% purity) were purchased from 
Praxair, Inc. and used as provided. For the gas-phase 
experiments, a sealable J. Young NMR tube was evac¬ 
uated to remove excess air, and filled with pure gas to 
15 PSIA. Liquids were purchased from Sigma-Aldrich. 
For liquid-phase experiments, samples were degassed and 
flame-sealed in NMR tubes. 


Nuclear Magnetic Resonance Methods 

Measurements were performed on a 14.1 T vertical 
bore Bruker AV 600 MHz NMR spectrometer equipped 
with a 5 mm broadband probe with a z-direction gradi¬ 
ent. The sample was placed in the center of the NMR 
magnet, and the pulse sequence (Fig. 1, main text) was 
applied. The receiver was placed in DQD mode (for¬ 
ward Fourier transform, quadrature detection of complex 
data). During acquisition, a linewidth-broadening mag¬ 
netic field gradient of amplitude g was applied. 
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Temperature Control 

Temperature is altered by the variable temperature 
(VT) unit, which controls the temperature of the sam¬ 
ple by heating or cooling the surrounding air. Tem¬ 
peratures above ambient temperatures (298 K < T < 
460 K) are achieved by the probe’s internal heater. 
Temperatures below ambient temperatures (200K <T 

< 298 K) are achieved with a heat-exchange coil, by 
pre-cooling the gas in a liquid nitrogen bath. Due 
to experimental limitations in the VT-system, we were 
unable to achieve temperatures lower than 200 K or 
above 460 K. In order to determine the real sample 
temperature, each VT-temperature value was calibrated 
against a neat methanol standard for temperatures 200 K 

< T < 300 K (Eq. 18) and a neat ethylene glycol stan¬ 
dard for temperatures 301 K < T < 460 K (Eq. 19): 

T[K] = 409.0 - 36.54(Ad) - 21.85(Ad)2 (18) 

and 

T[K] = 466.5 - 102.00(Ad), (19) 

in which (Ad) is the chemical shift difference between the 
two peaks of neat methanol or two peaks of neat ethylene 
glycol [23]. We note that at colder temperatures (j 290K), 
temperatures tended to fluctuate more than at ambient 
and above ambient conditions. 

SM: DATA ANALYSIS AND SAMPLE DATA 

Fast Fourier transformation (EFT) and subsequent 
analysis with the MATLAB Curve Fitting Toolbox was 
performed in MATLAB R2010b (The MathWorks, Inc.; 
Natick, MA). 


liquids, the lower bound of temperature was determined 
by the freezing point of the substance; all data above this 
temperature was fit to a linear regression. 

Figure 6(A) is the raw data from a methane linewidth 
vs. temperature experiment. The gradient strength is 
g = 0.15 G/cm, yielding a FOV of 78 cm. Figure 6(B) is 
the raw data from a dichloromethane linewidth vs. tem¬ 
perature experiment. The applied gradient strength was 
g = 0.3 G/cm, yielding a FOV of 94 cm. Linewidth 
for the gas-phase data follows a monotonic decrease as 
a function of temperature, whereas linewidth for liquid- 
phase data does not follow a clear trend. 

Linewidth vs. gradient data 

For linewidth vs. gradient data, the two regimes were 
determined based on the field of view (FOV), as calcu¬ 
lated according to the following equation: FOV=fs/^ng 
(fs, sampling rate; g, gradient amplitude). When the 
FOV becomes comparable to the size of the rf-sensitive 
region (~ 1 cm in our experiments), we are in the 
frequency-encoding (g^) regime; when FOV 3> 1 cm, we 
are in the non-frequency encoding, g^ regime. The coef¬ 
ficient of the g^ and g'^ slopes were determined by fitting 
data in the frequency-encoding regime and non-frequency 
encoding regimes respectively to a weighted quadratic 
function. 

Figure 7(A) is the raw data from a methane linewidth 
vs. gradient experiment. The data was acquired at 
T = 220 K. Figure 7(B) is the raw data from a 
dichloromethane linewidth vs. gradient experiment. The 
data was acquired at T = 300 K. Although the quadratic 
curves are drawn as a guide to the eye, the data illustrate 
that linewidth monotonically increases as a function of 
gradient strength. 


Linewidth vs. temperature data 


Determination of Error Bars 


For gases, linewidth vs. temperature data for each gas 
was fit to a linear regression. Only the high-temperature 
regime (i.e. above the Sutherland’s constant C) could be 
studied due to limitations in the VT control. In Ref. [21] 
the Sutherland’s constant was calculated for each gas ac¬ 
cording to its viscosity [21,24] and Sutherland’s model 
[21,25]: 


At each temperature point, five scans were separately 
acquired and averaged to reduce the effects of tempera¬ 
ture fluctuations over the time course of the experiment. 
For each substance studied, multiple experiments were 
performed on different days to eliminate the random fluc¬ 
tuations in field homogeneity day to day resulting from 
magnetic field drift. 


, /r\3/2 fr + c\ 
\f') (,r-f c j ’ 


( 20 ) 


SM: REFERENCES 


in which fj, is the viscosity of a gas at temperature T, g! 
is the viscosity of that same gas at another temperature 
T', and C is the Sutherland’s constant for that gas. The 
values of S were calculated to be 198 K, 237 K, and 292 K 
for methane, acetylene, and propylene respectively. For 


21. S. Ghapman and T.G. Cowling. The Mathematical 
Theory of Non-Uniform Gases: An account of the kinetic 
theory of viscosity, thermal conduction and diffusion in 
gases. Cambridge University Press, 1939. 

23. C. Ammann, P. Meier, and A.E. Merbach. J. Magn. 
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Chemistry and Physics, 95th Edition (2014-2015). 
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A) B) 
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FIG. 6. A) Raw linewidth vs. temperature data for methane gas. The applied gradient strength was g = 0.15 G/cm, with 
FOV = 78 cm. B) Raw linewidth vs. temperature data for liquid dichloromethane. The applied gradient strength was g = 0.3 
G/cm, with FOV = 94 cm. 
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FIG. 7. A) Raw linewidth vs. gradient data for methane gas, acquired at 220 K. B) Raw linewidth vs. gradient data for liquid 
dichloromethane, acquired at 300 K. The continuous curves drawn through the data points are guides to the eye, and are not 
actual fits. 



